Mon. Not. R. Astron. Soc. OOP. [TlfSTI C2009') Printed 30 November 2009 (MN WI^ style file v2.2) 



The nature of the Sloan Digital Sky Survey galaxies in 
various classes based on morphology, colour and spectral 
"q>^" features — III. Environments 

o : 
o ■ 

^ Joon Hyeop Lee^'^*, Myung Gyoon Lee^*, Changbom Park^*, Yun- Young Choi^* 

^ Korea Astronomy and Space Science Institute, Daejeon 305-348, Korea 

^Astronomy Program, Department of Physics and Astronomy, Seoul National University, Seoul 151-742, Korea 
Korea Institute for Advanced Study, Dongdaemun-gu, Seoul 103-722, Korea 

Astrophysical Research Centre for the Structure and Evolution of the Cosmos, Sejong University, Seoul 143-747, Korea 
Accepted 2009 November 23. Received 2009 November 23; in original form 2009 September 09 

ABSTRACT 

We present a study on the environments of the SDSS galaxies divided into fine classes 
based on their morphology, colour and spectral features. The SDSS galaxies are clas- 
sified into early-type and late-type; red and blue; passive, Hii, Seyfert and LINER, 
which returns a total of 16 fine classes of galaxies. We estimate the local number 
density, target-excluded local luminosity density, local colour, close pair fraction and 
the luminosity and colour of the brightest neighbour, which are compared between 
the fine classes comprehensively. The morphology-colour class of galaxies strongly de- 
pends on the local density, with the approximate order of high-density preference: red 
early-type galaxies (REGs) - red late-type galaxies (RLGs) - blue early-type galax- 
ies (BEGs) - blue late- type galaxies (BLGs). We find that high-density environments 
(like cluster environments) seem to suppress AGN activity. The pair fraction of Hii 
REGs does not show statistically significant difference from that of passive REGs, 
while the pair fraction of Hn BLGs is smaller than that of non-Hn BLGs. Hn BLGs 
show obvious double (red -I- blue) peaks in the distribution of the brightest neighbour 
colour, while red galaxies show a single red peak. The brightest neighbours of Seyfert 
BLGs tend to be blue, while those of LINER BLGs tend to be red, which implies 
that the difference between Seyfert and LINER may be related to the pair interaction. 
Other various environments of the fine classes are investigated, and their implication 
on galaxy evolution is discussed. 

Key words: galaxies: general - galaxies: evolution - galaxies: statistics - galaxies: 
elliptical and lenticular, cD - galaxies: spiral -- galaxies: active 



1 INTRODUCTION 

Galaxies in the Universe have a wide range of physical prop- 
erties such as morphology, spectra, mass, and so on. One of 
the major goals of the extragalactic astronomy is to under- 
stand what determines such properties of galaxies, and the 
first step toward the goal is to understand exactly how dif- 
ferent the various galaxies are. One efhcient way to com- 
prehend the nature of the diverse galaxies is to classify 
them using some criteria and to compare their statistical 
properties between the classes. In many previous studies, 
however, galaxies have been classified using simple schemes 
based on just one or two properties, which sometimes miss 
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detailed aspects of galaxy evolution. For example, several 
unu sual classes of galax i es such as blue early-type galax- 
ies (jPerreras et al. 20051: Lee et al. 20061) and passive sp i- 



ral galaxies (|Yamauchi fc Gotdl2004 llshigaki etall 120071 ') , 
or the fundamental parameter difference betwe en finely- 
divided classes in the same morphology class (|Lee et al.l 
l2007l : [Ch oi et al ]|2009 a'l are difficult t o investigate using sim- 
ple classification schemes. Recently, iLee et al] (|2008l . |2009| . 
hereafter Paper I and Paper II, respectively) presented stud- 
ies on the optical and multi-wavelength properties of the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) galaxies 
divided into fine classes based on their morphology, colour 
and spectral features. Those finely-classified SDSS galax- 
ies show distinguishable differences from each other in their 
photometric (both optical and multi- wavelength) and struc- 
tural properties, which are hardly found in the studies using 
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simple classification schemes. Based on the understanding of 
the basic properties of the various galaxy classes presented 
in Paper I and Paper II, now it is needed to approach the 
next question: what makes them have such diverse natures? 

For a long time, the environments of galaxies have 
been regarded as important factors determining the na- 
ture of galaxies. Since close galaxies have influence on 
each other directly or indirectly, the local density about a 

J alaxy is closely related to the propert ies of that galaxy 
Kuehn fc Rvdenll2005l : IPark et al.ll2007l ). Environments in 
which galaxies have many chances to interact with each 
other, induce morphology transfor mation by the tidal in- 
teraction or merging of galaxies |Spitzer fc Baadd Il95ll : 
iMerrittI [iQS^ : ICoziol fc Plauchu-Fravnl bOOTh . Such galaxy 
interactions con vert the properties of galaxies, by triggering 
star formation (IWoods et al.|[200^ or by supp ressing star 
formation with tidal gas stripping arziani e t al. 2003). 
Moreover, the complete merging of two or more galaxies 
is a classica l candidate of the early-type g alaxy formation 
mechanism (jToomrdl 19771 : ISearle fc Zinnl ll978). In addition 
to the direct interaction between galaxies, dense environ- 
ment s affect gala xies with their strong gravitational poten- 
tial (|Bekkilll999l ). Since large amount of hot gas is con- 
centr ated in high-den sity environments (e.g. galaxy clus- 
ters; [Ch urazov et al. 

I [1996), the morphology transforma- 
tion of galaxies ha ppens in such environments by the ram 
pressure stripping (Gunn fc GottI Il972l : iQuilis et al. I I2OOOI : 
iRoediger fc HensleIj|2005^ . 

Due to the close relationships between environments 
and galaxy properties, many studies have b een focused on 
the en vironmental effects in galaxy evolution. iBlanton et al.l 
l|2005h inspected the environmental dependence of optical 
properties of the SDSS galaxies, finding that the colour is 
the most sensitive parameter to the local density, and that 
the other parameters are not sensiti ve to environment s when 
the colour and luminosity are fixed. iKuehn fc RvdenI (|2005l ) 
found that the surface brightness profile type and axis ra- 
tio o f galaxies are fr e quent ly correlated with the local den- 
sity. iBernardi et al.l (120061 ) showed that the age, metallic- 
ity and a-element enhanceme nt of galaxies a re aff ected by 
the environments of galaxies. iMateus et al.l l|2007l ) argued 
that galax;y evolution is accelerated in dense environments, 
which took place mainly at high redshifts, based on the age - 
density relations wi th respect to galaxy luminosity or mass. 
[Rogers et al.l (|2007t ) investigated the star formation history 
of early-type galaxies, showing that a significant immber of 
galaxies in high-density environments undergo low but de- 
tectable recent star formation. Owers et al. (200^ inspected 
the environmental dependence of starburst galaxies, finding 
that the starburst in a galaxy is triggered when its nearby 
(within 20 kpc) neig hbour has similar luminosity or mass. 
IWestobv et al.l (|2007l ) presented that most AGN host galax- 
ies are closer to the red sequence than the blue sequence, but 
they pre fer lower density environments than red sequence 
galaxies. iPark et al.l (|2007f ) inspected the environmental de- 
pendence of physical properties of the SDSS galaxies, based 
on the local density estimated using an adaptive smooth- 
ing kernel, showing that variations of galaxy properties are 
almost entirely due to the environmental dependence of mor- 
phology and luminosity, and that other properties are almost 
independent of local density, when morphology and luminos- 
ity are fixed. Ivan den Bosch et al] (|2008l ) found that the role 



of satellite galaxies in the transformation from blue sequence 
galaxies to red sequence galaxies is important, in the sense 
that the satellite gala xies may quench the s tar formation 
of their host galaxies. IPark et all (|2008l ) and IPark fc Choi 
(2009) showed that the morphology of a galaxy depends 
sensitively on the distance and morphology of its nearest 
neighbour galaxy, and that galaxy properties in the general 
field do not directly depend much on the large-scale back- 
ground density. They showed that the apparent dependence 
of galaxy properties on the background density environment 
can be largely explained by an accumulated result of the in- 
teractions between the nearest neighbour galaxies. 

Since it is obvious that environments affect galaxy prop- 
erties significantly, environments must have played an im- 
portant role in determining the fine classes (i.e. morphol- 
ogy, colour, and spectral features) of galaxies as well. Thus, 
to understand the origin of the various fine galaxy classes, 
it is needed to compare their environments. This is the 
third in the series papers on the nature of galaxies in var- 
ious classes based on the morphology, colour and spectral 
features. In this paper, we investigate the environments of 
the finely-divided galaxy classes, using several environmen- 
tal parameters. The outline of this paper is as follows. Sec- 
tion 2 shows the data set we used, and §3 describes our 
classification scheme of the SDSS galaxies and the defini- 
tion of the environmental parameters. In §4, we present the 
results, and their implication on galaxy evolution are dis- 
cussed in §5. The conclusions in this study are summarised 
in §6. Throughout this paper, we adopt the cosmological 
parameters h = 0.7, SIa ~ 0.7, and Q,m ~ 0.3. 



2 DATA AND SELECTION CRITERIA 

We briefly summarise the datasets, classification scheme and 
sample selection in this paper. More details are described in 
Paper I and Paper II. 



2.1 Sloan Digital Sky Survey 

We use the SDSS Da ta Release 4 (DR4; 
lAdelman-McCarthv et al.ll2006l ) main galaxy spectroscopic 
sample. The SDSS is a photometric and spectroscopic sur- 
vey, mapping about one quarter of the whole sky, and the 
SDSS DR4 spectroscopic data cover 4783 deg^. The median 
FWHM of point sources in the r band images is 1.4", and the 
wavelength coverage in the spe ctroscopy is 3800 - 9 200A. 
We use both the SDSS pipeline (jStoughton et al.ll200 data 
and the Max-Pl anck-Institute for Astronomy catalogue 
(MPA catalogue; Kauffmann et all l2003al : iTremonti et al.l 
l2004l : lGallazzi et 81.112006^ Using the SDSS atlas images, we 
have measured the colour gradient, inverse concentration 
index, and the axis ratio of g alaxies corrected for the incli- 
nation and the seeing effects (|Choi et al.ll2007l ). Corrections 
of the SDSS magni tudes have been con ducted for the fore- 
ground extinct ion JSchlegel et al. Il998h . the redshift effect 
(K-correction; iBlanton et al.l l2003t ) and galaxy luminosit y 
evolution (evolutionary correction; iTegmark et al.l |2004| ). 
We corrected the observed magnitudes of about 360,000 
SDSS galaxies into the magnitudes at redshift z = 0.1, 
because the SDSS galaxies were observed most frequently 
there. Since the corrections are optional in this paper. 
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Figure 1. Morphology and colour class distribution on the BPT 
diagram for a small sample of the SDSS galaxies. Different sym- 
bols indicate different morphology-colour classes: REGs (filled cir- 
cle), BEGs (open circle), RLGs (filled star) and BLGs (open star). 
The line shows the boundary dividing Hll galaxies and AGN host 
galaxies in this paper. 



Figure 2. Absolute magnitude versus redshift for the SDSS 
galaxies. The solid-line box represents the sample volume selected 
in this paper. The dashed lines are the limits of the complete spec- 
troscopic sample. 

Table 1. Abbreviations of the 16 fine galaxy classes 
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we denote the magnitude with K-correction as '^'^^m and 
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Blue 


the magnitude with both K-correction and evolutionary 
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pRLG 


pBLG 
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Hii 
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hBEG 


/iRLG 
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sBEG 


sRLG 


sBLG 




LINER 


/REG 


/BEG 


/RLG 


/BLG 



2.2 Additional multi- wavelength datasets 

As a near- infrare d dataset, we use t he Two Micron All Sky 
Survey (2MASS: [jarrett et allbOOOD Extended Source Cat- 
alogue (XSC). The number of 2MASS extended sources 
matched with our SDSS sample is about 200,000, with 
false match probability of 0. 005 per cent. The InfraR ed As- 
tronomical Satellite fI RAS: iNeugebauer et al.1 ITsol ) Faint 
Source Catalogue (FSC: [Moshir et al.lll992l ) is used as a mid- 
and far-infrared dataset. About 7,000 IRAS sources were 
matched with our SDSS sample, and the false match proba- 
bility between SDSS and IRAS is about 0.071 per cent. We 
use the F aint Images of the R adio Sky at Twenty-centimetres 
(FIRST; iBecker et all Il995h data, and the number of the 
matched FIRST sources is about 14,000, with 0.45 per cent 
false mat ch probability. The Galaxy Evolution Explorer 
(GALEX: [Marziani et al] 120031 ) survey GR2/GR3 data are 
used Eis an ultraviolet dataset. About 85,000 GALEX objects 
were matched, with the false match probability about 1.46 
per cent. As an X-ray dataset, we use the Roentgen Satellite 
(ROSAT: lAschenbach et al.lll98ll ) all sky survey data. The 
information of the ROSAT sources matched with the SDSS 
objects provided by the SDSS is used, and the number of 
the ROSAT sources is about 2,200. 



2.3 Classifications and volume-limited sample 

In this paper, we classify the SDSS galaxies using three cri- 
teria: morphology, colour and spectral features. Here, we 
briefly introduce the galaxy classification scheme. More de- 
tails about classifications and possible biases are described 
in Paper I and Paper II. First, the SDSS galaxies were 
classified into early-type galaxies and late-type galaxies in 
the colou r - colour gradient - light-concentration parame- 
ter space (jPark fc Choill2005h . Second, we divided the SDSS 
galaxies into red galaxies and blue galaxies, based on the 
colou r distribution of early-type galaxies as a function of red- 
shift ()Lee et al.ll2006l ). Third, the SDSS galaxies were classi- 

Table 2. The number of galaxies in each fine class 





REG 


BEG 


RLG 


BLG 


Total 


Passive 


5109 


26 


446 


77 


5658 


Hii 


2479 


148 


1626 


9632 


13885 


Seyfert 


1477 


98 


1770 


1472 


4817 


LINER 


4829 


159 


3569 


1281 


9838 


Total 


13894 


431 


7411 


12462 


34198 
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fied into passive galaxies, Hii galaxies, Seyfert galaxies and 
low ionisation nuclear emission region (LINER) galaxies, 
based on their flux ratios bet ween several spectral lines (e.g. 
BPT diagram ; Baldwin et al.. 1981 : Kauffman n ct a l. 2003b; 
iKewlev et al.ll2006l ). Fig. [Tlshows the distributions of differ- 
ent morphology and colour classes on the BPT diagram. As 
mentioned previously, it is noted that the three kinds of clas- 
sification do not always agree with each other, showing the 
necessity of fine classification. Those three galaxy classifica- 
tions with different criteria return a total of 16 fine classes 
of galaxies: [early-type, late-type] x [red, blue] x [passive, 
Hii, Seyfert, LINER]. Hereafter, we use the following ab- 
breviations for the 16 galaxy classes: REG (red early-type 
galaxy), BEG (blue early- type galaxy), RLG (red late- type 
galaxy), BLG (blue late- type galaxy), and p- (passive), h- 
(Hii), s- (Seyfert), /- (LINER). The abbreviations for the 16 
fine galaxy classes are summarised in Table [T] For RLGs, 
we only use RLGs with their axis ratio larger than 0.6, be- 
cause many BLGs with large inclinations may be classified 
into RLGs due to their internal dust extinction l|Choi et al.l 
l2007l ). 

Since most properties of galaxies are known to be sen- 
sitive to their luminosity and redshift, the sample volume 
needs to be carefully controlled in statistical studies of galax- 
ies. Among the three volumes selected in Paper I and Paper 
II, we use the VI sample as shown in Fig. (2] which contains 
the largest number of galaxies and covers the brightest range 
of absolute magnitude among the three volumes: —22.92 ^ 
°-^^^Mp^t{r^ ^ -20.71 and 0.08 < z ^ 0.10. The number 
of galaxies in each fine class in this sample is listed in Table 
[21 and more details about the sample-volume selection are 
described in Paper I. 



3 ENVIRONMENTAL PARAMETERS 
3.1 Local density indicators 

There are two different ways in estimating the density of 
local environments, according to whether we define the lo- 
cal volume two-dimensionally (2D) or three-dimensionally 
(3D). In the 2D method, the projected local density is esti- 
mated within a fixed redshift interval centred on the target 
galaxy. On the other hand, in the 3D method, the local 
density is estimated using the 3D distance calculated by re- 
garding redshift difference purely as distance difference. It is 
known that 2D-based parameters are more accurate in high- 
density environments, while 3D-based parameters are mor e 
accurate in low-density environments (jCooper et ahlfioOSh . 
Since both methods have their own merits and demerits, 
we estimate both of them. In this paper, we adopt ±1000 
km s"""^ as the redshift interval for the 2D local density es- 
timation. More detailed comparison between the 2D- and 
3D-based environmental parameters are shown in i]All 

The most frequently used parameter as a local density 
indicator is local number density (pjv). When estimating 
Pn, we adopt a spline smoothing kernel (see ijA2|) . In ad- 
dition, we use the comoving distance to the n-th nearest 



^ The Petrosian absolute magnitude in the r band, after K- 
correction and evolutionary correction. 



neighbour as the kernel size, which means that the smooth- 
ing kernel size is adaptive to the local number density. The 
adaptive smoothing prevents oversmoothing in dense regions 
and preserves better the property of galaxy distributions 
than a smoothing with a fixed kernel size, as described by 
I Park et al.l (|2007l ) , who used the spline- kernel with adaptive 
size enclosing a fixed number of the nearest neighbour galax- 
ies. In this paper, the local number density pjv in Mpc""^ is 
defined as: 



Pjv = ^ fsp{dk/d„)/V 



(1) 



where n is an integer to define the kernel size dn as the 
comoving distance to the n-th nearest neighbour, and d^ is 
the comoving distance to the k-th nearest neighbour. fap(x) 
is the spline kernel function and V is the local volume in 
Mpc^. It is noted that the d„ and dfc are projected comoving 
distances in the 2D estimation. The local volume V is a 
cylinder in the 2D estimation, the height of which is 2000 
km s~^, while it is a sphere in the 3D estimation. From now 
on, we denote the 2D- and 3D-based local number densities 
as Pn and ^° pN, respectively. 

The local galaxy number density is the most frequently 
used local density indicator. However, the number density 
does not reflect the local environments perfectly, because 
the properties of the counted galaxies are not considered 
in that parameter, such as their luminosity and colour. For 
this reason, local luminosity density or loca l mass density are 
also used as the lo cal density indicators l|Park et al.ll2008l : 
iPark fc Choill2009l ). In this paper, we estimate local lumi- 
nosity density in the r band (pr), which is defined as: 



Pr = -2.5 log ( /sp(dfc/dn)xlO 



/V] (2) 



in the unit of absolute-mag Mpc~^, where °'^^M)(;(r) is the 
K-corrected r band Petrosian absolute magnitude of the k- 
th nearest neighbour. In the same manner, the local lumi- 
nosity density in the u band (p„) is defined, and the local 
colour is simply defined as pu — pr- The behaviours of the 
local density indicators in this section are described in i]A3l 
with more details. Throughout this paper, we adopt n = 5 
for all local density indicators, to focus on the small-scale 
environments of galaxies (see i]A3[) . We corrected the lo- 
cal luminosity densities by excluding the information of the 
target galaxy: target-excluded local luminosity density {pr') 
and target-excluded local colour [pj — pr')- The necessity 
of this correction is described in i]A4l 

Since our sample luminosity range is biased to the 
very bright side of galaxy luminosity distribution, there is a 
caveat that the estimated local luminosity density and local 
colour may not exactly reflect the local environments. That 
is, the local luminosity density depends on the fleld-to-field 
variation of galaxy luminosity function. Moreover, the local 
colour may be bluer if we use a sample with more fainter 
luminosity range, because g alaxy colours ten d to be redder 
with increasing luminosity (IChoi et al.ll2007h . This should 
be kept in mind throughout this paper. 

Fig. [3] presents how pjv and pr' are correlated with the 
clustercentric distance, based on the 2D and 3D estimation 
respectively. The clustercentric distance is defined as the 
projected comoving distance to the nearest known galaxy 
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Clustercentric distance [Mpc] Clustercentric distance [Mpc] 



Figure 3. (a) ^°pjv versus the clustercentric distance, (b) pr' versus the clustercentric distance, (c) ^°pjv versus the clustercentric 
distance, (d) pr' versus the clustercentric distance, jipf (= 2.469 X 10~^ Mpc~^) is the mean number density in our sample. The 
clustercentric distance is defined as the projected comoving distance to the nearest galaxy cluster within the ±1000 km s~^ redshift slice. 
In each panel, the solid line connects the median local density values at given clustercentric distance bins. 




2 2.5 3 1.5 2 2.5 
- p^' [abs.mag Mpc~^] 

Figure 4. The target-excluded local luminosity density versus 
target-excluded local colour relation, based on (a) the 2D esti- 
mation and (b) the 3D estimation. The lines connect the median 
Pu' — Pr' at given pr' bin. 

clustei0 within the ±1000 km redshift slice. At the clus- 
tercentric distance < 2 Mpc, ^°p,.' decreases rapidly as the 
clustercentric distance decreases, while the pr' variation in 

^ The coordinate list of the galaxy clusters was re- 
trieved from the NASA Extragalactic Database (NED; 



|http : // nedwww. ipac . caltech .edu/j l . 



the cluster regions is much smaller than the pr' variation. 
If we adopt 1-2 Mpc as the typical radius of galaxy clus- 
ters, the intermediate environments between galaxy clusters 
and fields approximately correspond to about ^°pr' ^ —17 

C°Pn/pn ~ 20) and 17.5 C^pn/pn ~ 60), where 

pjv {~ 2.469 X 10"'' Mpc""') is the mean number density in 
our sample. It is noted that the spectroscopic completeness 
in the central re gions of rich clus ters may be lower than 
that in the field jVoon et alT 2008V This effect causes the 
underestimation of the local density in very-high-density en- 
vironments, making the correlation between the local den- 
sity and the clustercentric distance in Fig.[3]less clear. Fig.|4] 
shows the p/ versus pj — pr' plot, an 'environmental colour- 
magnitude diagram'. Like the galaxy colour-magnitude di- 
agram, the environmental colour-magnitude diagram shows 
a sequence, in the sense that brighter pr' environments have 
redder pj ~ pr' ■ 



3.2 Close pair fraction 

IWoods fc Geltol l|2007l ) defined close pair galaxies as galax- 
ies with their recession velocity difference < 500 km s^^ 
and thei r projected distance < 50 kpc. Adopting the defi- 
nition of IWoods fc Gelled l|200'iD . we estimate the pair frac- 
tion in each fine class. We define the brightest neighbour 
as the brightest one among the neighbour galaxies satisfy- 
ing the close pair definition for a given target galaxy. The 
r-band absolute magnitude and u — r colour of the bright- 
est neighbour are denoted as ''■'^'^M(r)BrightcstNcighbour and 
"■^^{u - r)BrightostNoighbour respectively, the distribution of 
which will be compared between the fine classes, in the next 
section. It is noted that the close pairs in this paper do not 
necessarily mean that they are currently interacting pairs, 
which should be mor e strictly defined using the virial radii 
of the pair galaxies (|Park et al.l l2008l : IPark fc Choil |2009| : 
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Figure 5. (a) The dependence of the early-type class fractions on the 2D local number density, (b) The same as (a), but for the late-type 
classes, (c) The dependence of the early-type class fractions on the target-excluded 2D local luminosity density, (d) The same as (c), but 
for the late-type classes. Errorbars show the Poisson errors. The intermediate range between cluster and field environments is denoted 
as the shaded domain. 



IPark fc Hwaii3l2009l ). However, due to the practical diffi- 
culty in estimating the mass-to-light ratio of each fine class, 
we use this simple definition of close pairs, which selects pair 
galaxies that have relatively high chance of recent interac- 
tions. 

4 RESULTS 

4.1 Local density 

Fig. [5] and [6] show the dependence of each class fraction 
on the local number density and the (target-excluded) lo- 
cal luminosity density. The class fraction in Fig. [5] and [6] 
is defined as the number fraction of the galaxies in a fine 



class among the entire galaxies, at given local density bin. 
The REG fraction at high-density environments is much 
larger than tha t at low-den s ity en vironments. However, as 
reported by West obv et al.1 (j2007l ). it is found that AGN 
REGs (sREGs and iREGs) prefer slightly lower density en- 
vironments than non-AGN REGs (pREGs and /iREGs). In 
Fig. [5^, non-AGN REG fractions monotonically increase as 
the local density increases, while AGN REG fractions have 
a peak at ^° pn/pn ~ 60 C°Pr' ~ —18). At denser environ- 
ments than the peak density, the number fraction of AGN 
REGs does not increase any more or even decreases with 
large uncertainty as the local density increases. It is noted 
that the peak local density is very close to the intermediate 
local density between cluster and field environments (see 
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Figure 6. The same as Fig.fS] but for the 3D local density indicators. 



Fig. 131). This feature is hardly found in Fig. [S] but the 3D 

estimation is less reliable than the 2D estimation at high- 
density environments. Note that the Y-axis scales in Fig. [5] 
and Fig.|B]are very different between the panels (i.e. between 
the fine classes). A better view of the fine class composition 
at fixed local density is presented in Fig. [T] which compares 
the class fraction between the fine classes at given 2D local 
luminosity density bin. 

Non-AGN RLGs (pRLGs and /iRLGs) show a trend 
similar to non-AGN REGs, in the sense that their number 
fraction increases as the local density increases. The frac- 
tion of AGN RLGs (sRLGs and ZRLGs) shows a peak at 
intermediate-density environments and the number fraction 
of AGN RLGs decreases as the local density increases at 
high-density environments (^°pjv/pjv > 20 or ^°pr < —17). 
In the 3D local density, this feature is also found in the 
ZRLGs but not found in the sRLGs (Fig. |6j. Fig. [8] com- 



pares the number fraction of red non-AGN galaxies (pREG 
+ /iREG + pRLG + /iRLG), red AGN-host galaxies (sREG 
+ IREG + sRLG + /RLG), blue non-AGN galaxies (pBEG 
+ /iBEG + pBLG + /iBLG), and blue AGN-host galax- 
ies (sBEG + /BEG + sBLG + /BLG). In this figure, it is 
clearly shown that the number fraction of AGN-host galax- 
ies rapidly decreases as the local density increases at high- 
density environments (which approximately correspond to 
galaxy cluster environments), with significant difference be- 
tween red non-AGN galaxies and red AGN-host galaxies. 

In Fig. [5] and [6l no significantly different trend is found 
between pREGs and /iREGs, which implies that the local 
density is not directly responsible for the star formation in 
/iREGs. The fraction of AGN BEGs (sBEGs and ZBEGs) in- 
creases as the local density decreases, whereas the non-AGN 
BEG fraction seem to have a peak at intermediate-density 
environments with large uncertainty. The BLG fraction at 
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Figure 7. Comparison of tlie class fraction between the fine classes at given 2D local luminosity density bin. In each panel, the X-axis 
bins indicate the fine classes in order of: pREG, hREG, sREG, /REG, pBEG, hBEG, sBEG, /BEG, pRLG, /iRLG, sRLG, /RLG, pBLG, 
/iBLG, sBLG and /BLG. The errorbars show the Poisson errors. 



low-density environments is much larger than that at high- 
density environments. The pBLG fraction as a function of 
local density is less biased to low local density than the other 
BLGs, but the uncertainty is large. 

4.2 Median local luminosity density and the local 
colour 

Fig.|9]plots the median p/ and pu' — pr' of the 16 fine galaxy 
classes with their sampling errors. The sampling errors were 
estimated from 200-times random sampling tests for given 
sample size, which show how statistically significant the dif- 
ferences in the median pr' and pu' — pr' are between dif- 



ferent fine classes. The median p/ and pu' — p/ values, 
their sample inter-quartile ranges (SIQRs) and sampling er- 
rors are summarised in Table |3] and [l] In Fig. [O] even the 
fine classes in the same morphology- colour class have often 
significantly different environments. For example, pREGs 
and /iREGs have similar environments Pr' ~ —16.2 and 
"°(p„' - Pr') ~ 2.45), but AGN REGs (sREGs and ZREGs) 
have significantly different environments from those of non- 
AGN REGs by more than 5(j, in the sense that AGN REGs 
prefer less dense and bluer environments Pr' ^ —15.7 and 
"°(pu' - Pr') ~ 2.39) than pREGs and /iREGs. AGN BEGs 
also prefer less dense and bluer environments than /iBEGs 
with 2.3cr significance. 
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Figure 8. (a) The number fraction of red non-AGN galaxies 
(open circle), red AGN-host galaxies (filled circle), blue non-AGN 
galaxies (open triangle) and blue AGN-host galaxies (filled trian- 
gle) as a function of the 2D local number density, (b) The same 
as (a), but as a function of the 2D luminosity density. Errorbars 
show the Poisson errors. The intermediate pr' range between clus- 
ter and field environments is denoted as the shaded domain. 



Table 3. Median 2D environmental parameters for the fine classes 



Class 


2D / 
Pr 




'""{Pu'-P 


r') 


pREG 
/iREG 
sREG 
IREG 


-16.206 ±0.038 
-16.089 ±0.056 
-15.726 ±0.069 
-15.786 ±0.035 


(1.535) 
(1.463) 
(1.397) 
(1.374) 


2.459 ± 0.008 
2.452 ±0.011 
2.396 ± 0.015 
2.377 ± 0.009 


(0.282) 
(0.277) 
(0.292) 
(0.274) 


pBEG 
hBEG 
sBEG 
/BEG 


-15.347 ±0.507 
-15.443 ±0.217 
-15.020 ±0.259 
-15.062 ±0.206 


(1.330) 
(1.229) 
(1.463) 
(1.152) 


2.390 ± 0.111 
2.424 ± 0.044 
2.287 ±0.056 
2.267 ± 0.044 


(0.215) 
(0.245) 
(0.253) 
(0.296) 


pRLG 
hRLG 
sRLG 
/RLG 


-16.335 ±0.128 
-15.982 ±0.067 
-15.567 ±0.058 
-15.692 ±0.043 


(1.678) 
(1.425) 
(1.413) 
(1.355) 


2.498 ± 0.026 
2.446 ± 0.014 
2.378 ± 0.014 
2.371 ± 0.009 


(0.300) 
(0.275) 
(0.269) 
(0.276) 


pBLG 
hBLG 
sBLG 
/BLG 


-15.277 ±0.282 
-15.093 ±0.027 
-15.142 ±0.068 
-15.196 ±0.074 


(1.575) 
(1.291) 
(1.407) 
(1.257) 


2.257 ±0.062 
2.304 ± 0.006 
2.292 ± 0.016 
2.318 ± 0.017 


(0.335) 
(0.292) 
(0.297) 
(0.281) 



±values show the sampling errors, and (value)s show the SIQRs. 

For late-type galaxies, the 2D and 3D environmental pa- 
rameters show somewhat different trends from each other. 
For example, pRLGs prefer significantly higher local density 
(by 5.5a) and redder local colour (by 4.1a) than AGN RLGs 
in the 2D estimation, while the significance of such differ- 
ence is relatively small in the 3D estimation {2a and 3.1cr, 
respectively). The local density and the local colour of BLGs 
are almost the same between their different spectral classes 



Class 


3D. / 
Pr 




'""{Pu'-P 


r') 


pREG 
/iREG 
sREG 
/REG 


-16.923 ±0.035 
-16.986 ±0.049 
-16.624 ±0.065 
-16.697 ±0.035 


(1.349) 
(1.298) 
(1.289) 
(1.279) 


2.442 ± 0.008 
2.452 ± 0.013 
2.415 ±0.016 
2.378 ± 0.010 


(0.279) 
(0.275) 
(0.291) 
(0.280) 


hBEG 
sBEG 
/BEG 


— iO.UoO it U.40Z 

-16.398 ±0.211 
-16.117 ±0.250 
-16.073 ±0.214 


^1 . ZoZ J 

(1.296) 
(1.419) 
(1.290) 


z.4iy it u.iUo 
2.418 ±0.051 
2.320 ± 0.063 
2.344 ± 0.046 


U.zoo ) 
(0.228) 
(0.286) 
(0.302) 


pRLG 
/iRLG 
sRLG 
/RLG 


-16.721 ±0.102 
-16.523 ±0.060 
-16.490 ±0.053 
-16.579 ±0.041 


(1.283) 
(1.415) 
(1.394) 
(1.343) 


2.460 ± 0.027 
2.426 ± 0.016 
2.367 ±0.014 
2.386 ± 0.010 


(0.298) 
(0.279) 
(0.296) 
(0.277) 


pBLG 
/iBLG 
sBLG 
/BLG 


-16.059 ±0.307 
-15.997 ±0.024 
-16.124 ±0.062 
-16.243 ±0.064 


(1.351) 
(1.393) 
(1.453) 
(1.383) 


2.400 ± 0.068 
2.313 ± 0.007 
2.288 ± 0.016 
2.340 ±0.018 


(0.302) 
(0.300) 
(0.303) 
(0.284) 



±values show the sampling errors, and (value)s show the SIQRs. 

in the 2D estimation. However, they show distinguishable 
difference in the 3D estimation, in the sense that /BLGs 
prefer slightly higher local density and redder local colour 
than sBLGs (by 2.5a). As mentioned in tj3.1l the 2D param- 
eters show good performance at high-density environments, 
while the 3D parameters work well at low-density environ- 
ments. Thus, it seems that the environments of RLGs are 
more reliable in the 2D estimation, whereas those of BLGs 
are more reliable in the 3D estimation. 

pRLGs have environments similar to those of pREGs, 
but pBLGs prefer environments significantly (3.3cr) less 
dense than those of pREGs. The environments of pBEGs 
are intermediate between those of passive red galaxies and 
pBLGs, but their sampling errors are very large (< la). 
As mentioned in the previous section, /iREGs prefer high- 
density environments in spite of their current star formation. 
/iRLGs have the local density similar to those of /iREGs, 
and the environments of /iBEGs are less dense than those of 
/iRLGs by 2.4cr. It is noted that /iREGs, /iBEGs and /iRLGs 
have different median p/ but similar median — p, '- 
However, /iBLGs prefer significantly (> 2.7a) bluer envi- 
ronments than the other Hii classes. No large difference is 
found between Seyfert galaxies and LINER galaxies. Among 
AGN host galaxies, AGN REGs prefer the highest-density 
environments, and AGN RLGs have environments slightly 
(1.8cr) less dense than AGN REGs. AGN BEGs and AGN 
BLGs prefer similarly low-density environments. 



4.3 Close pair 

Fig. llOl shows the relationship between the pair fraction and 
the local number density of the fine classes. As expected, it 
is found that the fine class with higher median local number 
density tends to have larger pair fraction in Fig. 1101 For ex- 
ample, the pair fraction of /iREGs pN /pN ~ 0.9) is about 
2 per cent, while that of /iBLGs pn/Pn ~ 0.55) is about 
0.55 per cent. In addition to this overall trend, some fine 
classes seem to show distinguishable difference in their pair 
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Figure 9. The target-excluded local luminosity density versus target-excluded local colour relation of galaxies divided into the 16 fine 
classes. Each symbol shows the median pr' and pu' — pr' for each class and and the errorbars indicate their sampling errors estimated 
from 200-times random sampling, using (a) 2D local density indicators and (b) 3D local density indicators. The cross symbol shows the 
median pr' and pu' — pr' of the entire galaxies. 
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Figure 10. Close pair fraction versus (a) '^^ pN /pN and (b) pN / PN ■ The errorbars show the Poisson errors of the close pair fraction 
and the sampling error of the Pn/pn- The datapoints with their Poisson errors > 1 per cent are denoted using light grey, and those with 
0.5 ^ Poisson error < 1 per cent are denoted using dark grey. Zero per cent pair fraction is denoted using light grey. 



fraction at given median local number density. For example, 
the pair fraction of /iBLGs is half of that of AGN BLGs 
despite their similar local number density, which shows the 
possibility that the pair interaction is responsible for the 
AGN activity in BLGs. Some other class sets seem to show 
pair fraction difference with similar local number density, 
such as 'sRLGs - ZRLGs - /iRLGs', '/REGs - sREGs' and 
'/iREGs - 'pREGs', but their uncertainty is large. 

Fig. [11] presents what properties the brightest neigh- 



bours in each fine class have: absolute magnitude, colour and 
local number density. It is found that the brightest neigh- 
bours of red galaxies have a narrow range of their colour, 
while those of blue galaxies show scattered colour distribu- 
tions. This difference does not seem to correlate with the ab- 
solute magnitude or the local number density. Fig. [12] shows 
the "^'^^(it — r) colour distribution of the brightest neigh- 
bours in each fine class. Red galaxies show similar "'^^(u — 
7')BrightestNeighbour distributions rcgardlcss of their fine class: 
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Figure 11. (a) The colour-magnitude relation and (b) the colour- density relation of the brightest neighbours in each fine class. 




Figure 12. The colour distribution of the brightest neighbours in each fine class. 



a single peak at ''■^^(?i-r)BrightestNcighbour = 2.8 and a small 
blue tail. On the other hand, /iBLGs and sBLGs show quite 
different "'^^(ii — r)BrightostNoighbour distributions from those 
of the red galaxies: the °-^^{u — r)BrightostNoighbour distribu- 
tions of /iBLGs and sBLGs show double peak features, with 
the primary blue peak at °-^^{u — r)BrightcstNcighbour = 1-6 
and the secondary red peak at °-^^{u - r)BrightcstNcighbour = 
2.8. The secondary red peak colour of the ftBLG and sBLG 
brightest neighbours is similar to the single peak colour of 
the red galaxy brightest neighbours. It is noted that the blue 
peak of sBLGs is twice in size of their red peak, whereas 
ZBLGs show a single red peak. Table[5]lists the pair fraction 



and median "'^^(ii — r) colour of the brightest neighbours in 
each class. 



4.4 Multi-wavelength detection and environments 

In Paper II, we presented several multi-wavelength proper- 
ties of the SDSS galaxies divided into fine classes. In the 
multi-wavelength study, the detection trends of the fine 
classes vary largely according to the survey wavelength. 
For example, the most/least efficiently detected classes are 
ZRLGs/pBEGs in the 2MASS, while those in the FIRST 
are sBLGs/pRLGs. Such variation is mainly due to the dif- 
ferences in the internal properties between the fine classes. 
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Figure 13. The dependence of the SDSS galaxy detection fractions in the 2MASS (circle), IRAS (pentagon), FIRST (rectangle), 
GALEX (triangle) and ROSAT (star) on (a) the 2D local number density, (b) the 2D local luminosity density, (c) the 3D local number 
density and (d) the 3D local luminosity density. The line shows the ROSAT detection fraction of the optically passive galaxies. The 
intermediate pj^ and pr' ranges between cluster and field environments are denoted as the shaded domains. The errorbars show the 
Poisson errors. 



but some galaxies show multi-wavelength properties that 
are not easy to understand just by their intrinsic charac- 
teristics, such as optically passive but X-ray bright galaxies. 
Those unusual properties are sometimes related to their ex- 
ternal conditions, because some environments (e.g. galaxy 
clusters) themselves have strong emissi on in some wave- 
length such as X-ray l|Bourdin et aL|[2004[ ). Thus, in this sec- 
tion, we briefly investigate the environmental dependence of 
multi-wavelength detection fractions for our sample galax- 
ies, which helps us to understand the multi-wavelength be- 
haviours of the fine classes better. 

Fig. [13] shows the dependence of the detection frac- 
tions in the multi- wavelength surveys on the local luminosity 



density. The 2MASS detection fraction increases as the lo- 
cal density increases, whereas the detection fractions in the 
IRAS and GALEX decrease with increasing local density. 
Considering the detection trends of the fine galaxy classes 
in the multi- wavelength bands (Paper II), Fig. 1131 confirms 
the previous finding that galaxies in high-density environ- 
ments tend to have red colours, old ages and early-type 
mor phology, compared to those in low-density environment s 
(|Blanton et al.ll2005l : iBernardi et"aLll2006l : iPark et al.lboOTi ') . 
The FIRST detection fraction shows quite different trends 
between the 2D and 3D local density parameters. The 2D lo- 
cal density of the FIRST detection fraction increases almost 
monotonically as the local density increases, while its 3D 
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Table 5. Close pair information for the fine classes 



Class 


pair fra.ction^'^^ 


0.1KC„ _ 

^ ' BrightcstNcighbour 


pREG 


1.51 lb 0.17% 


2.86 lb 0.15 


/iREG 


1.98 lb 0.29% 


2.82 lb 0.09 


a-rv.-Cjvj 


1 09 -1- n 9fi% 
-L.UZ ^ U.ZU/0 






i.oy it U.io/o 


Z.oO it U.l ( 


pBEG 


0% 






1.00 m u.yD /o 


9 Q9 -1- n f^n 


sBEG 


1.02 ± 1.03% 


2.38 ±0.00 


ZBEG 


2.52 ±1.27% 


2.21 ±0.37 


pRLG 


2.01 ± 0.68% 


2.90 ±0.05 


/iRLG 


0.86 ±0.23% 


2.76 ±0.21 


sRLG 


1.30 ± 0.27% 


2.73 ±0.22 


m,LG 


0.90 ±0.16% 


2.81 ±0.18 


pBLG 


2.60 ±1.86% 


3.58 ±0.96 


/iBLG 


0.54 ± 0.08% 


2.07 ±0.57 


sBLG 


1.29 ± 0.30% 


1.89 ±0.54 


ZBLG 


1.25 ±0.31% 


2.75 ±0.29 



(a) Close pair fraction. ±value shows the Poisson error. 

(b) Median O j.) colour of the brightest neighbours. ±value 
shows the SIQR. 



local density shows a peak at ^° pN / pN ~ 60 {^^ pr 



^18). 



Considering that the 2D local density is more reliable than 
the 3D local density at high-density environments, it seems 
like that the galaxies at higher-density environments have 
higher probability to be radio sources. 

The ROSAT detection fraction shows a very strong de- 
pendence on local luminosity density. In Fig. 113b .. the largest 
(at ^°pn/pn ~ 1700) and the smallest (at ™ pnIpn ~ 5) 
ROSAT detection fractions are different by 2.3 dex. This re- 
sult indicates that the bias of the ROSAT-detected galaxies 
to red colour (Paper II) is closely related to the environments 
of galaxies, because red galaxies prefer high-density environ- 
ments as shown in Fig. [5] In Fig. 1131 the ROSAT detection 
fraction of the passive galaxies (solid line) is very similar 
to that of the entire galaxies (open circles), showing that 
the passive galaxies in high-density environments have high 
probability to be detected in the X-ray. Since galaxy clusters 
are typically strong X-ray emitters with a l arge amount of 
intra c luster hot gas jMcMillan et al.| Il989l : iRhee &: Latourl 
Il99ll : iBourdin et ai]|2004 ). some ROSAT-detected passive 
galaxies in high-density environments may be affected by 
such cluster environments. However, some passive galaxies 
in low-density environments are also detected in the X-ray. 
The number fraction of the ROSAT-detected passive galax- 
ies at pr' > —17 among the entire ROSAT-detected pas- 
sive galaxies is about 30 per cent (16 of 54), which are classi- 
fied as X-ray Bright Optically Normal Ga laxies (XBONGs; 
lYuan fc Naravanil2004l : ICivano et al.ll2007l ). 



5 DISCUSSION 

5.1 Red Early- type Galaxies 

The environments of REGs are previously well studied, 
in the sense that galaxies at high -density envi r onments 
tend to have early-type morphology (|Oemleij| 19741 : iDresslei] 



Il980l : iPostman fc Gelleilll984l ) . Recently, iPark et all (|2007l ) 

showed that the number fraction of early-type galaxies 
monotonically increases as the local density increases at 
given luminosity. In this paper, pREGs and /iREGs show 
good agreement with those previous results, displaying 
strong and monotonic dependence of their number frac- 
tion on their local density. However, unlike those non-AGN 
REGs, the number fraction of AGN REGs (sREGs and 
ZREGs) is not a monotonic function of the local density. 
Instead, they have a peak at ^^pn/pn ~ 60 pr' ~ —18), 
which is very close to the intermediate local density be- 
tween cluster and field environments shown in Fig. O This 
difference between non-AGN and AGN REGs results in 
the difference in their median local luminosity density and 
local colour (Fig. [U]): AGN REGs prefer less dense and 
bluer environments than pREGs and /iREGs by more than 
5(7 significance. This result is closely related to the fact 
that AGN host galax;ies prefer the green valley domain in 
the colour- magnitude diagram, whi ch is intermediate be- 
tween the red and blue sequences ^Silverman et al.l l2008l : 
IChoi et aLll2009bl:lHickox et al.ll200 9''l . because galaxy colour 
and local density are tightly correlated (Park ct al. 2007). 

If we suppose that galaxies at very-high-density envi- 
ronments are mostly cluster galaxies (supported by Fig. 
cluster environments seem to suppress AGN activity, reduc- 
ing the number fraction of AGN RE Gs at very-high-density 
environments. iPark fc Hwand (|2009| ) discussed several pos- 
sible mechanisms causing a sudden variation of the galaxy 
morphology in cluster environments, such as galaxy-galaxy 
tidal interaction, galaxy-cluster potential interaction and 
galaxy-galaxy hydrodynamic interaction. Those mechanisms 
probably affect also the AGN activity in cluster environ- 
ments: the mechanisms such as galaxy-galaxy interactions 
or harassment may deprive cluster galaxies of cold gas that 
can be the fuel of AGNs, resulting in a sudden decrease of 
AGN host galaxies in very-high-density environments. 

During the study of the early-type star formation his- 
tory, [Rogers^eT^l] (|2007h found that a significant number of 
galaxies in high-density environments undergo low but de- 
tectable recent star formation. /iREGs in our sample may 
correspond to those objects, which show no significant dif- 
ference in the local density distribution from pREGs. Such 
similarity of environments between pREGs and /iREGs im- 
plies that the local density is not directly responsible for 
the star formation in /iREGs. Since satellite or close pair 
galaxies play an impo rtant role in determining the properties 
of the target galaxy l|van den Bosch et al ] |2008l : |Park' et al.l 
l2008l : iPark fc Choil l20oi r it is a possible scenario that the 
star formation in /iREGs is triggered by the interaction with 
their close neighbour galaxies. In Fig. 1101 it is found that 
/iREGs tend to have higher chance of pairs (1.98 ± 0.29 per 
cent) than pREGs (1.51 ± 0.17 per cent). However, this dif- 
ference is statistically not meaningful (1.4cr significance) in 
our sample. It is an open question why the AGN-suppressing 
mechanisms (discussed in the previous paragraph) seem not 
to suppress the star formation in /iREGs. 

5.2 Blue Early- type Galaxies 

Since the number fraction of BEGs out of the entire galaxies 
is very small (about 1.3 per cent in our sample), there are 
not many previous statistical studies on their environments. 
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In IPark et all (|2007l ). it has been shown that the number 
fraction of BEGs out of early-type gal axies decreases as 
the lo cal density increases. More recentlvjKannapp an et aP 
found that many blue-sequence E/SO galaxies reside 
in low-to-intermediate-density environments. Similar results 
are shown in this paper that the BEG fraction at low-density 
environments is higher than that at high-density environ- 
ments. However, there seems to be a difference between non- 
AGN BEGs and AGN BEGs, in the sense that AGN BEG 
(sBEG and ZBEG) fraction increases monotonically as the 
local density decreases, whereas the non-AGN BEG fraction 
has a intermediate-density peak (Fig. [S] and O. This is sim- 
ilar to the relationship between non-AGN REGs and AGN 
REGs. However, due to the very small sample size of BEGs, 
the uncertainty is so large that the statistical reliability of 
the results for BEGs is not high (less than 2.3(j significance 
from Fig. O. 

In Paper I and Paper II, BEGs show the photomet- 
ric and structural properties implying that they are early- 
type galaxies in the forming phase. Since it is known 
that high-density en vironments accelerate galaxy evolution 
l|Mateus et ahllioOTl ) , the preference for low-density environ- 
ments of BEGs seem to be natural, explaining why those 
early-type galaxies (BEGs) are forming so lately compared 
to REGs. In detail, however, the environmental discrepancy 
between non-AGN BEGs and AGN BEGs is not easy to 
explain, because AGN-host galaxies are thought to be in 
the intermediate phase fro m star-forming galaxies to passive 
galaxies (AGN feed back; lAntonuccio-I)elogu fc SUM l2008l : 
iRaflertv et allbOOSl ). 

5.3 Red Late- type Galaxies 

RLGs are mostly regarded as bulge-dominated late-type 
galaxies (Paper I and II). RLGs are not frequently men- 
tioned in the previous studies, but several studies pro- 
vided some hints on their environments. That is, it 
is known that galaxies have greater bulge-to-disc ratio 
(or high light-conc entration) with increasing local density 
i Giuricin et al.lll995. : Hashimoto fc Oemler.l999i). As prob- 
able bulge-dominated galaxies, RLGs agree with those pre- 
vious findings. Non-AGN RLGs (pRLGs and /iRLGs) show 
a trend similar to non-AGN REGs, in the sense that their 
number fraction increases as the local density increases. On 
the other hand, the AGN RLG (sRLG and ZRLG) fraction 
shows a peak at intermediate-density environments, which 
is a trend similar to that of AGN REGs. In the same con- 
text as AGN REGs (SJSHI, it seems that the AGN activity 
in RLGs are suppressed at very-high-density environments. 

The unanswered problem in i]5.1l (whv very-high-density 
environments do not suppress the star formation?) is less 
serious in RLGs, because the number fraction of ftRLGs 
slightly decreases at the highest-density environment with 
large uncertainty (Fig. O 1.2a significance). However, this 
decrease is statistically not meaningful and the star forma- 
tion in /iRLGs seems to be clearly less suppressed than the 
AGN activity in AGN RLGs at high-density environments. 
To understand this problem better, further studies on the 
Hii galaxies and AGN host galaxies in cluster environments 
are needed in the future. 

In Paper I and II, pRLGs have properties very similar to 
those of REGs. Interestingly, in this paper, pRLGs have en- 



vironments similar to those of pREGs. This result supports 
the idea that pRLGs are intermediate objects between early- 
type galaxies and late-type galaxies as suggested in Paper I 
and II. 

5.4 Blue Late-type Galaxies 

BLGs mostly consist of spiral or irregular galaxies, which 
have the opposite trend to REGs in their environments. 
That is, galaxies with late- t ype morphology tend to be 
at low -donsity environments (Oeinie? 'l974l: iDressle^ Il980l : 
iPostm an fc GcUcr 1984; Park ot al. 2007). In this paper, 
consistent with those previous findings, the number fraction 
of BLGs increases almost monotonically as the local density 
decreases. Exceptionally, pBLGs seem to prefer intermediate 
density environments, but the sample size is too small. 

In Fig. 1101 /iREGs have a similar pair fraction to 
pREGs. On the contrary, /iBLGs have a smaller pair fraction 
than the other BLGs (significant up to 2.4(7) despite their 
similar local number density. This result indicates that the 
close pair interaction is less important for the star formation 
of BLGs compared to that of REGs, at least in our sample 
luminosity range (-22.92 <o ike Mpet{r) < -20.71). On 
the other hand, the close pair interaction may be responsi- 
ble for the AGN activity in BLGs. 

In Fig. 1121 the different spectral classes of BLGs show 
impressive differences in their close pair colour distribution. 
/iBLGs have both a blue peak (~ 1.6) and a red peak 
(~ 2.8) in their ° "^^(m - r)BrightostNoighbour distribution, and 
the two peaks have similar size. sBLGs also show the dou- 
ble peak feature, but the red peak is just half in size com- 
pared to the blue peak. On the other hand, iBLGs have 
only a red peak, which is similar to that of red galaxies. 
The difference in the brightest neighbour colour between 
sBLGs and iBLGs also affect the median local colours of 
them, in the sense that sBLGs prefer slightly bluer local 
colour than /BLGs (Fig. [SJa). It is noted that the single red 
''■^^(ii — r)BrightcstNcighbour peak of ZB LGs is unusual, consid- 
ering the conformity of pair galaxies l|Weinmann et al"]|2006l : 
Par k et al.,,200^ ). These results show a possibility that the 
supply of fuel gas from outer sources is closely related to 
the AGN type: Seyfert or LINER. In other words, at least 
for BLGs, it seems like that Seyfert galaxies tend to have 
blue close pairs that may supply rich fuel gas for the target 
galaxy, unlike LINER galaxies. 

However, this conclusion is very cautious for two rea- 
sons. First, the uncertainty in the °'^^(u — r)BrightostNoighbour 
distribution of BLGs is large, because the number of BLGs 
hosting close pair galaxies is not large enough. Second, such 
a trend (the difference in the "'^^(k — r)Brightes5tNeighbour dis- 
tribution between Seyferts and LINERs) is not found in the 
other morphology-colour classes. This may be due to the 
limited range of the sample galaxy luminosity, but it is not 
assured in this paper. Nevertheless, it is not easy to neglect 
the difference between sBLGs and /BLGs shown in Fig. [TT] 
and 1121 either. This problem needs more inspection in the 
future, using a larger sample of close pairs. 

5.5 Comparison in a given spectral class 

pRLGs have environments similar to those of pREGs, but 
pBLGs prefer environments significantly less dense and bluer 



Environments of SDSS galaxy classes 15 



than those of pREGs {3.3a significance). The local environ- 
ments of pBEGs are intermediate between those of passive 
red galaxies and pBLGs, but their sampling errors are large 
(less than la significance). The preference for low-density 
environments of pBLGs is not surprising, because many of 
them may not be genuinely passive (Paper II). Except for 
pBLGs, the other passive galaxies prefer relatively dense 
and red environments as expected. All ROSAT-detected 
optically-passive galaxies are passive red galaxies (Paper II) 
and about 70 per cent of them reside in high-density envi- 
ronments like galaxy clusters. Since a galaxy cluster itself 
is a strong source of X-ray, some of those galaxies may be 
contaminated by the galaxy cluster X-ray emission. 

In spite of their current star formation, /iREGs pre- 
fer high-density environments. It is interesting that /iREGs, 
ftBEGs and /iRLGs have different pr' distributions but sim- 
ilar pu' — pr distributions. On the other hand, /iBLGs pre- 
fer significantly (more than 2.1a significance) bluer environ- 
ments than the other Hii classes. Since the luminosity of 
our sample galaxies is limited to the bright end, this does 
not necessarily mean that Hii galaxies except /iBLGs do not 
have blue neighbours. If /iREGs, /iBEGs and ftRLGs have 
some blue neighbours, the neighbours may be fainter than 
the lower luminosity limit (°-^^^Mpet(r) = -20.71) of our 
sample. The pair fraction of ^REGs is similar to that of 
pREGs, while the pair fraction of /iBLGs is smaller than 
that of non-Hii BLGs (significant up to 2.4(7). This indi- 
cates that the role of close pairs is less important for /iBLGs 
than for /iREGs. 



6 CONCLUSIONS 

Here, we summarise the new findings in this paper and their 
implication on galaxy evolution. 

(1) The morphology-colour class of galaxies strongly de- 
pends on the local density. The approximate order of high- 
density preference is REGs - RLGs - BEGs - BLGs. 

(2) In low-density environments, the fraction of AGN host 
galaxies increases as the local density increases, but it de- 
creases in high-density environments that approximately 
correspond to cluster environments. This indicates that 
high-density environments (like cluster environments) sup- 
press AGN activity. 

(3) The local luminosity density distribution of /iREGs is 
almost the same as that of pREGs. It is an open question 
why the AGN-suppressing mechanisms at high-density en- 
vironments do not suppress the star formation in /iREGs. 

(4) The non-AGN BEGs prefer intermediate-density envi- 
ronments, while the AGN BEGs prefer low-density environ- 
ments like BLGs, although with low significance. In the con- 
text of early-type formation, such preference for low-density 
environments of BEGs seems to refiect the later formation 
of early-type galaxies at lower-density environments. 

(5) pRLGs that have properties very similar to those of 
REGs (Paper I and II) have environments similar to those 
of pREGs, which confirms the idea that pRLGs are inter- 
mediate objects between early-type galaxies and late-type 
galaxies as suggested in Paper I and II. 



(6) The pair fraction of /iREGs does not show statistically 
significant difference from that of pREGs (1.4(t significance), 
while the pair fraction of /iBLGs is smaller than that of non- 
Hii BLGs (significant up to 2.4(t), indicating that the close 
pair interaction is less important for the star formation of 
BLGs compared to that of REGs. 

(7) Red galaxies show a single red peak in the distribu- 
tion of the brightest neighbour colour, whereas /iBLGs and 
sBLGs show double (red -I- blue) peaks. The brightest neigh- 
bours of sBLGs tend to be blue, while those of ZBLGs tend 
to be red, which shows the possible relationship between 
AGN types and close pair interactions. 

(8) Some ROSAT-detected passive galaxies in high-density 
environments may be affected by galaxy cluster environ- 
ments. However, about 28 per cent (15 of 54) of the ROSAT- 
detected passive galaxies reside in low-density environments, 
which are classified as X-ray Bright Optically Normal Galax- 
ies (XBONGs). 
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APPENDIX A: LOCAL DENSITY ESTIMATION 
Al 2D and 3D Estimation 

Since the location of a galaxy in the Universe is described 
in a three-axis frame at given time, wo should estimate 
the number density of galaxies in the 3D space, in princi- 
ple. Practically, however, the distance estimation for distant 
galaxies has very large uncertainty compared to the determi- 
nation of the 2D coordinates on the plane perpendicular to 
the lino-of-sight, due to the intrinsic limit of the distance es- 
timation. The distances to galaxies are typically measured 
by means of determining the recession velocities of galax- 
ies caused by the cosmological expansion. However, it has a 
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Figure Al. The difference between the 2D and 3D local density 
estimates as a function of the 2D local density: (a) local number 
density and (b) local luminosity density. The different line type 
indicates the different redshift interval in estimating the 2D local 
density: Av = ±500 km (long-dashed line), ±1000 km 
(solid line) and ±2000 km s^^ (short-dashed line). 



problem that the recession velocity of a galaxy is given by 
both the cosmological expansion and the peculiar motion of 
that galaxy. 

The velocity dispersions of typical rich galaxy cluster s 
range from 500 km s"^ to 2000 km s~^ (|Teague et al.lll990l ). 
Suppose that a galaxy cluster at z = 0.1 has velocity dis- 
persion of 1000 km s~^. In that cluster, 32 per cent of its 
member galaxies have relative velocities larger than 1000 
km s~^, compared to the dynamical centre of the cluster, 
which results in redshift difference larger than 0.003. If we 
estimate the distance to the member galaxies with ignoring 
the peculiar velocity effect, the 0.003 redshift diflerence is 
interpreted as the distance difference of 13 Mpc. Since the 
typical diameter scale of galaxy clusters is at most a few 
Mpc, this causes significant underestimation of the 3D local 
number density in rich galaxy clusters. On the other hand, 
since such a peculiar velocity effect is very small in the low 
density fields, the 3D local density is more accurate than the 
2D local density in low-density environments. 

Fig- ED compares between the 2D and 3D local density 
indicators, showing a clear trend that the 3D local density 
is smaller than the 2D local density at high-density environ- 
ments, while the 3D local density is larger than the 2D local 
density at low-density environments. The difference between 
the 2D and 3D densities tends to increase as the redshift in- 
terval in estimating the 2D local density increases. 



A2 Smoothing Kernel 

We adopted a spline smoothing kernel when estimating the 
local density, because it is centrally weighted unlike the uni- 
form kernel and because it has a finite size unlike the Gaus- 
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Figure A2. Smoothing kernels for local density estimation as a 
function oi d/d„: uniform (short-dashed) and spline (solid), d is 
the distance from the target object, and d„ is the kernel size. 



sian kernel. The basic functional f orm o f the spline kernel is 
described well in [cib ezon et al.l (|2008l ). In this paper, the 
spline kernel is defined as: 



fsp{d/d„) = A ■ 



2(1 - d/d„f 




^ d/d„ s; 0.5 
0.5 < d/d„ < 1 
d/dn > 1 



where d is the distance from the target object and d„ is the 
kernel size. A is a constant, which has the value of 40/7 and 
32/3 in the 2D and 3D local density estimation, respectively. 
Fig. IA2I shows the shape of the spline kernel used in this 
paper. 

A3 Local Density Indicators 

As a local density indicator, the local galaxy number den- 
sity is quite useful but has some weakness at the same time. 
That is, in the local number density estimation, both five 
faint blue galaxies and five bright red galaxies are commonly 
counted as five galaxies, which is missing the luminosity and 
colour information of the counted galaxies. The luminosity 
dispersion problem can be minimised by limitin g the lumi- 
nosit y of the counted galaxies to a narrow range (|Park et al.l 
l2007h . However, another way to reflect both the luminosity 
and colour information of the galaxies in the local environ- 
ments is to estimate the local luminosity density and the 
local colour. It should be reminded that the local luminosity 
densi t y is not the same as the local mass density (|Park et al.l 
l2008l : IPark fc Choill2009l ). because galaxies with the same 
mass may have quite different luminosity according to their 
mass-to-light ratio. Thus, both the local number density and 
the local luminosity density need to be compared for more 
reliable environmental analysis. 

Fig. lA3b , shows the distribution of the local galaxy num- 
ber density (pjv) divided by the mean number density in our 
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Figure A3, (a) The distribution of the 2D local number density (^°pjv) divided by the mean number density (pjv) for the SDSS 
galaxies in the volume-limited sample. The solid line, short-dashed line and long-dashed line show the results using ds, dio and 6.20, 
respectively, (b) The 2D local luminosity density (^°Pr) distribution in the r band, (c) Local colour (^°(pu — Pr)) distribution, (d) The 
local number den sity - local lumin osity density relation, using d^. The line shows the least-squares fit using the ordinary least-squares 
bisector method l|lsobe et al.lll99(t . 



sample, pN = 2.469 x 10"^ Mpc"^. Fig. lA3b and Fig. lASl 
are the same as Fig. lA3b . but for the local luminosity density 
(pr) and the local colour (p„ — pr), respectively. When the 
larger n value of the kernel size dn is adopted, the distribu- 
tions of pn/pn, Pr and p„ — pr become shallower, showing 
that a large kernel tends to smooth out the variations of 
small-scale environments. Since the main interest of this pa- 
per is the local environments that affect directly the proper- 
ties of galaxies rather than the large-scale structure, we focus 
on the small-scale environments by adopting the ds kernel. 
Fig. IA3t l and IA4t i display the local number density versus 



the local luminosity density in the 2D and 3D estimation, 
showing tight correlations with the ordinary least-squares 
bisector (|lsobe et al.lll990l ') fit: 

= -2.635 ±0.002 log(™p]v/pjv) - 14.731 (A2) 

^°Pr = -2.643 ± 0.003 log(^°p]v/pjv) - 14.654 (A3) 

with p,.-direction rms of 0.33 and 0.35 absolute-mag Mpc~'^, 
respectively. That is, the local luminosity density reflects 
well the local number density overall, but there are typically 
Pr dispersions of ±0.33 (±0.35) mag Mpc~^ at flxed pjv in 
the 2D (3D) estimation. 
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Figure A4. The same as Fig. IA3I but for tiie 3D local density indicators. 



A4 Corrected Environmental Parameters 

The local luminosity density and the local colour provide ad- 
ditional information that the local number density does not 
show. However, there is one issue we should consider: pr and 
Pu — Pr include the information of the target galaxy itself. In 
the local number density, this is not a significant problem, 
because the target galaxy itself is always counted as simply 
one galaxy, which affects the zeropoint of the environmen- 
tal parameter but does not cause any systematic distortion. 
On the other hand, this may cause a significant distortion 
for the local luminosity density and the local colour, because 
target galajcies have their own luminosity and colour ranging 
widely. 

Suppose a target galaxy and its several neighbour galax- 



ies and we calculate the local luminosity density and the lo- 
cal colour for the target galaxy. After that, we replace the 
target galaxy with a new one with quite different luminosity 
and colour. Then, the new local luminosity density and lo- 
cal colour may be different from the old values. However, the 
local environment has not changed actually, but it is only 
the target galaxy that changed. Therefore, it is reasonable 
to exclude the information of the target galaxy itself when 
estimating the local luminosity density and the local colour, 
to investigate the environmental effects on the tar get galaxy. 
The local number or ma s s den sity used bv |Park"e t al. (2003, 
l2008h and lPark fc Choil l|2009l ') excluded the contribution of 
the target galaxy itself. 

Fig. IA5I and I A6I compare the original environmental pa- 
rameters {pr and Pu — pr) and the target-excluded environ- 
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Figure A5. (a) The comparison between tlie original 2D local luminosity density and the target-excluded 2D local luminosity density. 
The line shows the ordinary least-squares bisector fit. (b) The dependence of the original 2D local luminosity density on the target 
luminosity, (c) The dependence of the target-excluded 2D local luminosity density on the target luminosity, (d) The dependence of the 
original 2D local colour on the target luminosity, (e) The dependence of the target-excluded 2D local colour on the target luminosity, (f ) 
The dependence of the original 2D local colour on the target colour, (g) The dependence of the target-excluded 2D local colour on the 
target colour. 



mental parameters {pr and pj — pr')- pr is systematically 
fainter than pr (because only four galaxies are used in the 
luminosity summation, rather than five), but still correlated 
with Pr (Fig.lXSh andlMkV In Fig.Hll3,d,f andH6}3,d,f, it 
seems that p,- and pu — pr strongly depend on the luminosity 
and colour of the target galaxy. However, such dependency 
becomes very weak when the pr and p„ — pr are replaced by 
Pr and pj ^ pr' , the target-excluded environmental param- 
eters, as shown in Fig. lASb .e.g and lA6fc .e.g. After excluding 
the information of the target galaxy, the environmental pa- 



rameters have correlations with the target luminosity and 
colour, as follows: 



= 0.484 X " '^MpetCr) - 5.171 



- Pr') = 0.101 X "■''^(u - r) + 2.091 



(A4) 
(A5) 



= 0.480 X Wlp^tir) - 6.002 (A6) 
'°(pn' - Pr') = 0.077 X "■''^(u -r) + 2.151, (A7) 
which show relatively weak dependence on the local envi- 
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Figure A6. The same as Fig. IA5I but for tiie 3D local density indicators. 



ronments, compared to the target-included local density in- 
dicators. 



